


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1953-05 


One-dimensional analysis of steady-flow 
air-water mixtures In pipes 


Booth, Merson. 


Massachusetts Institute of Technology 


http://ndl.handle.net/10945/24783 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
sa Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 
KNOX appointed — and published — scholarly author. 





LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 


http://www.nps.edu/library Monterey, California USA 93943 





ONE-DIMENSIONAL ANALYSIS 
Oli Uke eee we 
MIXTURES IN PIPES 
MERSON BOOTH 
1953 





Library 
U.S. Naval Pos,-raduate Schoct 
Monterey, Califorma 

















ONE-DIMENSIONAL ANALYSIS OF STEADY-FLOW AIR-\ATE 


MIXTTRES IN PIPES 


Long 
VERSON BOOTH, LIEUTENANT, ". 3. NAVY 
P.S., Ue. Se NAVAL ACADEMY 


OLAS 


SUPMITTED IN PARTIAL FULFILLMSNT 
OF THE REQUIREMENTS FOR THE 
TEGREE OF NAVAL ENGINEER 
at the 


MASSACITISETTS INSTITITE OF TECHNOLOGY 


May, 1953 








- 


pia Eas 


ONE-DIMENSIONAL ANALYSIS OF STRANY={FLOW AIR~WATER 
MIXTURES IN PIPES 


bbs 
Merson Booth, Lieutenant, ". S. Navy 


Submitted to the Nepertment of Naval Architecture and “erine Fneineering On 
May 25, 1953, in partial fulfillment of the requirements for the degree of 
Naval Engineer. 


ABSTRACT 


The object of this thesis is to investigate a mumber of simple flow 
exarples of air-water mixtures ané to present the results in « simple quick 
reference form which can be valuable for engineering use in design and under- 
standing the phenomena of flow of two-component mixtures. The problems 
analyzed sare: 


l. The one-dimensional trajectory of droplets accelerated in «a 
gan stream 


PPR eet Rok MME Ds 1 Boao pak) Mehl hel Mas) Mob ae) tC REC Dab ge 
tion process 


3. The variation of stream properties due to vell friction with 
water present in the etrean,. 


The procedure used was entirely analytical. The results of the first 
problem were calculated, tatulated, and plotted in dimensionless form for 
various water and air velocities, droplet diameters, distance moved down the 
duct, and water and sir properties. The plots should be valuable for engi- 
neering use in determining the one-dimensional trajectory of spherical perti- 
cles in 4 ges stream. An iteration procedure using the plots can be used vhere 
air velocity, droplet diameter, or stream properties vary during the acceleration. 


The results of the second ani third problem were calouleted, tabuleted, 
and plotted in dimensionless form for various water to air velocity ratios, 
water to air masa rates, air Mach mumbers, and stream property ratios. The 
results should be valuable for engineering use in all regions except near 
choking conditions in the stream. Near this condition the air velocity increases 
too rapidly for the droplets to accelerate with the air. Since this ves one 
of the assumptions in the analysis the plots are not correct in this rogion. 


Wiel wet hs eM MRA Mb ee Mr) kU eee gab) wee Oe Oe 1) oles tae 8 
term accounting for the effect of momentum exchange of droplets with the duct 
wall. The exact nature of this effect is not know, but some correlation with 


experimental data is given. 
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Z2__INTRODUCTION 


The object of this thesis is to investigate a mmber of simple flow 
examples of air-water mixtures. The problems to be analysed earet 
(1) the one-dimensional trajectory of droplets 
accelerated in a gas stream, 
CM lsh Mel ll at ce mL 
acceleration process, and 
(3) the variations of stream properties due to 
wall friction with water present in the stream. 

This investigation is undertaken as part of the Aerothermopresscr 
Project at Massachusetts Institute of Technology. The Aerothermopressor 
is a device to increase the stagnation pressure of a gas stream by reducing 
the stagnation temperature through evaporation of water injected into the 
flow. 

In the work on the Aerothermopressor, there has been no quick refer- 
ence guide to aid in design or in understanding the phenomena of flow of 
air—vweter mixtures. The purpose of this thesis is to analyse these separate 
effects and to present the results in a simple quick reference form which 


can be valuable for engineering use in problems concerning flow of air~ 
water mixtures. 
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bie telet- Toho oa| 
The example studied is the injection of one droplet at a given initial 


velocity into e moving gas stream. The problem is to determine the subse- 
quent droplet velocities as eae function of distance moved downstream and the 
properties of the gas and droplet. 
The assumptions made are: 
ODMR ct me kh ee leh a te 
CRC tak eee oot Va tal bet fT 
diameter during the acceleration, and 
(3) the table of drag coefficient va. Reynolds 
number as given in reference (1) and repeated 
pO Mey. 3- Bee Melee | 
The details of this analysis are given in the appendix. It was necessery to 
perform a graphical integration of a function of Reynolds mmber. The results 
Sao) Me SMT) ME helen beet EM bh MBE Mb a4 5a): MO OPE 4s0 > nt 
gration was done using the trapozoidal rule. The results of both these inte- 


grations are given in Table II. 


i TRF 
Figures I and II are non-dimensional plots of droplet velocity vs. dis- 


tance moved downstream as a function of gas velocity, droplet diameter, end 
the properties of the gas anc droplet. This same information is given in 
Table III. Figure I is for droplet velocity less than gas velocity; Figure 
II is for droplet velocity greater than gas velocity. The plots in Figures I 
and II are normalized to specific initial droplet velocities. For droplet 
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Velosity less than ges velooity, the initial droplet velocity is assumed zero. 
Por droplet velocity greater than gas velocity, the initial droplet velocity 
is assumed twice the air velocity. for an initial droplet velocity other 
than the normalized conditions, it is necessary to subtract the distance 

read at the initial condition from that of the final condition. This is 
shown in more detail in the appendix. 

Figures I and II are used in the following manner: From given initial 
droplet velocity, gas velesity, droplet diameter, and gas and droplet proper- 
ties, calculate me rT 6 “ay « Enter the plot with ie quantities to deter- 
mine the initial oe ye intersection of ordinate a Neotel Molt aa 
corresponding to a ocean lanai Tee o a art 
eR “— At eh 4 curve to the intersection of that curve and the 
OTT Ba ie oo al - PE from the initial point. Read off the velue of 
xf a eres oo the distance from the initial point with the required a or 
ae v. ACM Pet Bay. sk Gm : 

Figures ae and IV are plots of distances required to reach specific 
acW As) | a ue = .5,.1 as a function of droplet diameter, air Mach mmber, and 
ek) bey See ee assuming initial droplet velocity is sero. Figure III 
is for air stagnation temperature of Vi ae. stagnation pressure of 1 atmos- 
phere. Figure IV is for air stagnation temperature of 1500°R and stagna- 
tion pressure of 1 atmosphere. 

If the gas velocity, droplet diameter, or the etream and droplet proper~ 
ties are not constant, the solution can be approximeted by assuming the varie- 
tion to oecur in a stepwise fashion and apply the plots successively over 
each step. This can be done in the case of evaporation in which e curve of 


droplet diameter vs. droplet velocity or distance may be assumed, or in the 





case of a converging or diverging section in which, because of known area 
changes, 4 curve of gas velocity ve. distance can be assumed. 


Conclusion 
The results obtained are given in dimensionless form in es simple 4 


manner as possible. The plots should be valuable fer engineering use in 
determining one-dimensional trajectory of spherical particles in a gee strean. 
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Zit. DISCONTINUITY ANALYSIS 
Procedure 


The example studied is the flow of an air-water mixture past two sections 
1 and 2 in the flow. The problem is to determine the stream properties at 
section 2 as a function of the water velocity at 2, the water rate, and the 
stream properties at section 1. 
The aseumptions made are: 
(1) Constant area flow 
CRE beste REM cle hein temperature 
CR rts a te aie 
(4) No change in temperature of droplets 
(5) Perfect gas relation holds for air. 
The details of this analysis are given in the appendix. The general method 
of performing the analysis is given in reference (2). 


— T is Vv 
SO i One e eee e e 
various water mas eo MB! Cols mB elise 5 i “The star condition is a4 aaa 
condition corresponding to oY = 0 and M = 1 4nd does not depend on the water 
be O71. BAD) Me ba ayer in Tables IV through VIII. The guiding prin- 
> GSES Mt bk.) MbAsl MLE) mee BR As) he =o, is constant in going between tables. 
Figures V through X show the soo of droplet st bat 2 1 
properties assuming that at section l, a = 0, and at section 2, v. = 1 for 


numbers on eAch curve indi- 
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various initial Mach mombers and water rates. The 
i] 

cate the water rate to air rate ratio, — , ami the letter S on some branches 
ry 

of the curves indicate that a normal shock in the flow is necessary to reach 


bAet-t | Me be ae) 





ae 


Figures V through X are used in the following mannert Knowing the 
RO oa) od Ok ee Le RE Lh aA, Ee Go) ee hot a 
the required plot with the M, and Read on the ordinate the value of the 
pa ROE Re Mel -L sha Mba) a ae hr sd ASE 2 to that at section 1. 


Conclusion 

The plots given are easy to use and are correct within the assumptions. 
However, near choking conditions the air velocity increases very rapidly. It 
is erroneous to assume that the droplets will be able to acoelernte with the 
stream in the absence of very large drag coefficients. Therefore, near 
choking conditions the plots will not give a correct answer to the variation 
of stream properties due to droplet acceleration. However, the values of el 
where choking of the flow at section 2 is indicated is approximately correct. 
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Mesa 
The example studied is the flow of an eir-water mixture in a pipe with 
wall friction and droplets striking the wall. The problem is, knowing the 
initial stream properties and friction factor, to determine the strean proper— 
ties at any section downstream. 
The assumptions made are: 
(1) Constant area flow 
(2) Adiabatio—no change in stagnation temperature 
(3) Mo evaporation of water droplets 
(4) Mo change in temperature of droplets 
(5) Perfect gas relation holds for air 
CPR b et te the eC ers Phar TTT] ot ae 
(7) oe 
(f) In each differential section an amount of water leaves 
the stream and hits the wall,and an equal amount of 
water leaves the wall and in picked Ph oR os Mme ME bo te 
(9) That this amount of water is proportional to the 
length of section 
The details of this analysis are given in the apperdix. The general method 
eee Mek ee oe ee ene) 


Results 


ap it leo Soe eh ak 
Calculated for various water er aes numbers. The star condition is 
® normalised condition corresponding to choking of the flow. These values, 
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to slide rule seowrncy, ore tabulated in Tables XI] through XV and plotted in 
Mieures XI through XVII. These fiqures show the effect of friction and 
droplet momentum exchange with the walls on the stream properties. The mzbers 


w 
on each curve injicate the water rate to air rete ratio, m9 ° 


Ee 
The symbol ¥ in the term G + r| Lex incorporates as a pseudo-fric- 


tional term the effect of momentum exchange of droplets with the wall. This 
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RCE Ra Ca a Ce tte te ee ma Meo ee eee 
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and the droplet diameter and perhaps other atream properties. The nature 

of this function is not known. UWowever, using the data of reference (3) 

for sand together with the low Mach number analysis described in the appen- 

dix, it appears that K= .4 = for 1" diameter pipe and sand diameter of 

200 » and 450 pe. It appears that the constant of proportionality is slightly 

smaller for sraller sani diameters. The effect of duct diameter is no known. 
Figures XI through XVII are used in the following manner: (Knowing the 

initial stream properties enter Ficure XI with the initial Mach mumber, iL 


Se Ase Ba au Then knowing Fc 


the length of duct from the initial section 1 to any other seotion 2, calculate 


Gs Cr ee Cnet Gz J Nei, Gs oon Gz 4 aa 
Enter Figure XI with Eg ‘ ae 4 cee and s Send int off Moe Then entering 
OC tka tC ke ee re eee os 
the normalized property at sections 1 and 2. ie this form the ratio of the 
property at 2 to that at 1. Knowing the value of the property at 1, the v-lue 
Coe ae toe eo BARE Re ok 
Mgure XVII is used to determine the Mach number efter 6 normal shock 
knowing the Mach murrber and -, just, before the shock and assuming me = 1 on 
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TCO aL ee te Ct ak ees Creat. = ste ures Cr 
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point. Follow across horisontally with constant a until reaching the same 


Ww 
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Conolusion 

Provided more information can be obtained concerning the pseuco-fric- 
tional term, K, the plots given can be valuable for engineering use. This 
could well form the basis of a thesia to determine whether this simple model 
of the flow is correct and to correlate data in the field to two-component 
weet) 

In addition, very near choking conditions this model of the flow is not 
correct since the air velocity increases very rapidly in a very short length. 
It is erroneous to assume that the dronlets can accelerate with the stream 
in the absence of very large drag coefficients. Therefore, near choking 
comitiona the plots will not give a correct answer to the variation of 


stream properties due to friction. 
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A=CROSS SECTIONAL AREA 
C= SPEED OF SOUND 
Co= DRAG COEFFICIENT 
Cp= SPECIFIC HEAT AT CONSTANT PRESSURE 
D= HYDRAULIC DIAMETER OF DUCT 
dur= DROPLET DIAMETER 
f=FRICTION COEFFICIENT OF DUCT 
F= IMPULSE FUNCTION (PAt FAV?) 
h = SPECIFIC ENTHALPY 
K= RATIO OF SPECIFIC HEATS (Ce/cv) 
Lmax= MAXIMUM LENGTH OF DUCT FOR CHOKING FLOW 
M= MACH NUMBER 
m= MASS OF DROPLET 
P= STATIC PRESSURE 
fo" P,= ISENTROPIC STAGNATION PRESSURE 
R = REYNOLDS NUMBER 
T = TEMPERATURE (ABSOLUTE) 


os i ad aa STAGNATION TEMPERATURE (ABSOLUTE) 
= TIM 


V = VELOCITY 
w= MASS RATE OF FLOW 
X = DISTANCE ALONG DUCT 
P= MASS DENSITY 
Y = KINEMATIC VISCOSITY 
( *= REFERS TO A NORMALIZED CHOKING CONDITION 
()a= REFERS To AIR 
( )w= REFERS TO WATER 
( )o = REFERS TO ISENTROPIC STAGIVATION CONDIT/ON 
()1 = REFERS TO SECTION | 
()2= REFERS TO SECTION 2 
()rR= REFERS TO CONDITION RELATIVE TO AIR 





DROPLET TRAJECTORY-DETAILS OF ANALYSIS 
INDIVIDVAL DROPLETS JNJECTED INTO MOVING AIR STREAM 


LET: Va = CONSTANT 
DeFiNe: Va Va na (For Va > Vw) 
Vee Vue-Va (FoR VaX< Vu) 
Pa Ps Vi 


At i look 
FORCE EQUATION ON DROPLET 
mdVe _ _cy fe A Vr 
dt 2 
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INTERMS OF REAL DISTANCES 
X,2 = Valiz + Xayz (3) 
NOTE: IF Va > Vw, XR /5 IW NEGATIVE DIRECTION 
IN EQUATION (4) SIGN OF LAST TERM /s: 


+ IF Va < Vw 
— JF Va>p> Vu 
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NORMALIZING EQUATION (4) SO THAT INITIALLY, 
AT POINT], Va= Va ; SO THAT: 


Vw =O FOR Va > Pn 
Vur = 2Veq FOR EN 4 Vur 


LET THIS NORMAL/ZED POSITION BE POINT O, AND 
ANY OTHER POINT DOWNSTREAM BE POINT X. 
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DISCONTINUITY ANALYSIS 


EFFECT OF DROPLET ACCELERATION ON AiR STREAM 


ASSUMPTIONS: 
. ABIABATIC- NO CHANGE IN STAGNATION TEMPERATURE 
2. NO CVAPORATION OF WATER DROPLETS 
3.NO CHANGE JIN TEMPERATURE OF WATER DROPLETS 
4.PERFECT GAS RELATION WOLDS FOR AIR 
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DISCONTINUITY ANALYSIS - CON 'T 
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MODIFIED FANNO LINE ANAYS/S 
FRICTIONAL, ADIABATIC, CONSTANT-AREA. 


ASSUMPTIONS : 

|. ADIABATIC - WO CHANGE JN STAGNATION TEMPERATURE 

2.N0 EVAPORATION OF WATER DROPLETS 

3.NO CHANGE IN TEMPERATURE OF WATER DROPLETS 

4. PERFECT GAS RELATION HOLDS FOR AIR 

Me Lila tM LV a a 

6. FOR DROPLET ON THE WALL; Yt-0 

7-IN EACH OIFFERENTIAL SECTION, LENGTH dx, AN 
AMOUNT OF WATER, dW, , LEAVES THE STREAM 
AND HITS THE WALL. /N ADDITION AN AMOUNT OF 
WATER , dure , LEAVES THE WALL AND [5 PICKED 
uP BY THE STREAM. 


CONTINUITY: 
PV = CONSTANT 
ree + dV = 0 (1) 
‘4 Yj 
EQUATION OF STATE: 
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DEFINITION OF M: 
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MODIFIED FANNO LINE ANALYSIS - CON'T 
PIAA 


Wi Na GP ALLY 
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IN STFADY STATE WITH WATER OW THE WALLS, 
ASSUME THE FILM THICKNESS !I5S CONSTANT WITH 


TIME AT ANY SECTION AND THE FILM VELOCITY 
IS NEGLIGIBLE. 
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MODIFIED FANNO LINE ANALYSI5- CON T 
MOMENTUM: 








MOMENTUM 
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ASSUME THAT IN ANY SPECIFIC CASE THE AMOUNT OF 
WATER PICKED UP BY THE STREAM FROM THE WALLS 

IS) PROPORTIONAL TO THE LENGTH. 


PE. A 


RESULTS OF SIMILTANEOUS SOLUTION OF EQUATIONS 
(1) THROUGH (5) ARE PRESENTED IN TABLE X 





MODIFIED FANNO LINE ANALYSIS- CON'T 
TABLE X 
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INTEGRATION OF EQUATIONS : 
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FRE SULTS OF INTEGRATION AND INTEGRAL RELATIONS 
ARE GIVEN /N TABLE XT 











MODIFIED FANNO LINE ANALYSIS - CON’T 
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FORMULAS IN INTEGRAL FORM 
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MODIFIED FANNO LINE ANALYSI5 - CON’) 


VALUES OF 4E +k] max; 3, a a ‘ae 


2 ey eee 
ARE GIVEN FOR VALVES OF M FROM O702 ; 
iLL es el OL = 27 


XY AND PLOTTED IN FIGURES XI THROUGH XVI. 





MODIFIED FANNO LINE -LOW VELOCITY ANALYSIS — 
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Frictional, Adiabatic, Constant=-Area Flow 


Perfect Gas, K = 1.4, _ ee 
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